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N-body simulations of dark matter halos show that the density is cusped near the center of the halo. The 
density profile behaves as r~ 1 in the inner parts, where 7 ~ 1 for the NFW model and 7 ~ 1.5 for the 
Moore's model, but in the outer parts, both models agree with each other in the asymptotic behavior of 
the density profile. The simulations also show the information about anisotropy parameter f3(r) of velocity 
distribution. f3 »* in the inner part and (3 Ri 0.5 (radially anisotropic) in the outer part of the halo. We 
provide some distribution functions F(E,L) with the constant anisotropy parameter f3 for the two spherical 
models of dark matter halos: a new generalized NFW model and a generalized Moore model. There are two 
parameters a and e for those two generalized models to determine the asymptotic behavior of the density 
profile. In this paper, we concentrate on the situation of (3(r) = 1/2 from the viewpoint of the simulation. 
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1. Introduction 

The dark matter halo can be considered as the collisionless self-gravitational system which is de- 
scribed by the Vlasov equation and the Poisson equation. Apart from observations, there are mainly 
two approaches to the study of dark matter halos: (1) numerical simulations and (2) analytic or 
semi-analytic method. 

Nowadays, N-body simulations become more and more important in the study of dark matter 
halos. Those simulations provide density profile and other properties of dark matter halos such as 
anisotropy of velocity distribution. A 'universal' density profile of dark matter halo was proposed by 



Navarro, Frenk, & White (hereafter NFwj22l. The NFW profile is shown as p oc l/((r/r s )(l+r/r s ) 2 ), 
where r s is a characteristic radius. We can find directly from the NFW profile that p oc r" 1 near the 
center and p oc r~ 3 in the outer parts. Some simulations' results exhibit different inner logarithmic 
slopes from NFW's but deviate not so much from the latter's in the outer parts of the halo. In the 
inner parts of the halo, Moore 's profikP behaves as p oc r 15 . Jing and Suto's profile^ behaves as 
p oc r~ 7 in the inner part, where the density logarithmic slopes 7 ~ 1.1 ~ 1.5 with different merger 
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histories and different total mass of the halos. Some other simulations concentrate on the outer slope 
of the dark matter density. For example, Avila-Reese et al. found that outer slopes 7 of some halos 
are larger than NFW's outer slopes(7 = 3^P. What's more, Hansen and Moore^^ found that the 
density logarithmic slopes 7 is correlated with the velocity anisotropy which is parameterized as an 
anisotropy parameter (3 and they provided the formula (3 ss 1 — 1.15(1 — 7/6). So (3 ~ in the inner 
part as 7 w 1 and (3 « 0.5 in the outer part as 7 ~ 3. 

On the other hand, many authors try to construct the dynamical models of the stellar sys- 
tem and the dark matter halo analytically or semi-analytically. It is important to construct the 
dynamical models which have physical meaning, as for instance, these analytic or semi-analytic dy- 
namical models can be used to generate the initial conditions of N-body simulations^! A stellar 
system or a dark matter halo is described by the distribution functions (DFs) F(x,v). Eddington 
(1916) showed how to determine the DF of a spherical symmetric stellar system with the isotropic 
velocity distribution^!, but it is difficult to calculate the distribution func tion in the anisotropic 
cases. A pioneer work on the anisotropic cases is called King-Michic modePEl that comes from 
an approximate steady state solution of the Fokker-Planck equation. In the past few years, there 
has been great progress in the anisotropic cases. Dejonghe to ok a la rge step on finding anisotropic 
distribution function with using the augmented density p(ib, r jHESI. Some authors constructed the 
anisotropic models wit h constan t anisotropy parameter^-H For another kind of models called the 
Osipkov-Merritt mod elpHll2l, the velocity dispersion tensor becomes isotropic near the center 
and becomes completely radial anisotropic in the outer part of the halo. CuddeforcP^ constructed 
the modified Osipkov-Merritt model which allows the velocity dispersion tensor to be an arbitrary 
anisotropy in the inner part of the halo but still completely radial anisotropy in t he outer part 
and the composite Osipkov-Merritt model was constructed by Ciotti & Pellegrini^. Furthermore, 
Baes & Hese^ constructed a dynamical model with seven parameters very recently to get a flexible 
anisotropy profile. Some authors built a model^ to erase the cusp in the center of the halo. 

To construct the dynamical models (or DFs) with flexible anisotropy beha vior, s ome special 
potential-density pairs have been consid ered. F or example, the Plummer 

model™, 

anisotropic 

Veltmann modepS, the Hernquist mo and the 7- 

mode |26 | 27|28 | 29 | 

Some authors study 

this matter in other ways. S ome ass um ed simple DFs at first and then solved the the potential 
and the density profile er^M Widrow^ provided the Osipkov-Merritt DFs and DEDs (differential 
energy distributions) for the NFW profile by the semi-analytic method. 

Recently, Evans and have provided the DFs with constant anisotropy of the dark matter 

for two types of the density profiles. One is the generalized NFW profiles, p oc l/(r(a + J") 6-1 ), 
and the other is the Gamma model, p cx l/(r 7 (a + r) 4-7 ). Both in these two models, asymptotic 
behaviors of density profiles are controlled by only one parameter. In this paper, we provide the DFs 
with constant anisotropy of the dark matter for more generalized density profiles in a semi-analytical 
way. Our generalized density profiles, with two parameters a and e, can cover many realistic profiles 
which come from simulations. 

In Section 2 we review the basic knowledge needed for this paper and the main ideas of DF with 
constant anisotropy. Section 3 provides DFs of two models for dark matter halos: a new generalized 
NFW model and a generalized Moore model. We make the discussion and conclusion in Section 4. 
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2. Basic Properties 
2.1. General formulae 

For the spherically symmetric stellar system with the isotropic velocity field, a mass distribution 
function F(E) describes this system very well. The mass density can be obtained from the distribu- 
tion function^, 



p(iP) = An / F{E)y/2{1> - E)dE , (1) 
Jo 

where the binding energy E is defined as 

E = ^{v)-\v 2 r -\v\ , (2) 
and ip(r) is the relative gravitational potential which can be obtained from the Poisson equation: 

r z ar \ dr J 
v r is the radial velocity and vt is the tangential velocity: 



v T = sj v e + v l ■ ( 4 ) 
EddingtorJSI provided the inversion formula of the Eq. |T]): 

f(f\- 1 n f E dpW d ± w 

f(e) -^ De Jo # vm^T)' ( ) 

where De denotes the differentiation with respect to E. 

The anisotropy parameter 21 mentioned in the Introduction is defined as: 

where a\ and of are the tangential and radial velocity dispersion. If j3 < 0, Oj, > 2af (tangentially 

anisotropic). If < j3 < 1, erf. < 2a^ (radially anisotropic). Else If j3 — 0, <r\ — 2 of (isotropic case). 

Specially, [3 — > — oo means that every particle is in a circular orbit and = 1 indicate that every 

particle is in a radial orbit. 

The anisotropic case is dif ferent from the isotropic one as we should consider the modulus of the 

angular momentum vectoif^l L = tvt and therefore now the distribution function of dark matter 

depends on two variables E and L. Then, the mass density can be obtained from F(E,L) by the 
rroi 

double integration^ 

r) = 2tt dE f^~ E) ,. F(E ^ , dvl . (7) 
Jo Jo \J2(^j -E)-v% 



However, the inversion formula of the above equation is much more difficult to obtain than that in the 
isotropic case. Many works have been done for this problem just as mentioned in the Introduction. 
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2.2. Distribution function with constant anisotropy 

There is a simple and widely used ansatz for F(E, L) as below-^^- 4 ! 18 ^^^ 4 ^ 

F(E,L)=L-^f(E) . (8) 

The simulations indicate that the anisotropy parameter 0{r) varies radially , but for the DF above, 
/3(r) is constant. Although this simple ansatz is very attractive in the calculation of inversion formula 
of Eq. the assumption of constant anisotropy parameter is not so ideal. Anyway, we can still use 
this valuable assumption which pave the way for the further work. 
From Eq. (|7|) and Eq. (j8]) we can get p from F(E, L), 

P&> r) = r-^ [% - E)V*-e m d E . (9) 



2/^(3/2-/3) 7 

The function f(E) can be given from the inversion of the above equation 

2^(2tt)- 3 / 2 _d_ f E d*P d n h 

r(i-A)r(i-/3) dEj (e - v^) A 

where h(ip) = r 213 p is expressed as a function of ip, n — [(3/2 — 0)\ and A = 3/ 2 — — n is the 
integer floor and the fractional part of 3/2 — 0. The DED now can be expressed as^UHl 

where te is defined through ^(je) = E. It is clear that te is the largest radius reachable by a 
particle with binding energy E. 

The expressions of DFs and DEDs are relatively simple if is a half integer (i.e., = 
3/2, 1/2, —1/2...). Specifically, in the case of = 1/2, which is more suitable from the viewpoint of 
the simulation, the expressions of DF and DED further reduce tcB3 



F{E:L) = 9±El d_M^ 

\ ' > 2tt 2 L ' dE Eyy h 

p + r(dp/dr) 

9{te) = 



(12) 



In this paper, we select 0=1/2 for our models just as Ref. EH by the same token. 

3. Dark Matter Halo Models 
3.1. New generalized NFW profiles 

Let us consider a family of density profiles with parameters a and e 

1 

(r/r s ) a (l +r/r s ) s 

We set the characteristic radius r s = 1, the total mass M tot = 1 and the gravitational constant 
G = 1 here, and then the density profiles reduce to 

P = C—~ -r^r , (14) 

r a (l + r) 3 + £ ~ Q 



P = C , , . . r -r= X . (13) 
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_ r(3 - a + e) 
° ~ 47rr(3 - a)r(e) ' ( &) 



With two parameters a < 3 and e > we can freely select the asymptotic behavior of the density 
profile both in the inner and outer part of the halo. For the cases of < e <C 1, the profile closes to 
the NFW profiles when a — > 1 and is similar to the Moore or Jing & Suto's profile if a = 1.1 ~ 1.5. 
If e = 1, the profile just reduces to the 7 model. 

We can calculate the relative potential ip from the Poisson equation. 



1&(r) =C[co- 4nr- a T(~a) (A^r) - 2A 2 (r) + A 3 (r))] 

At(r) = afitl-a + e, -a, 2- a, -r] (16) 
A 2 (r)= 2 F 1 [2-a + e, -a, 2 - a, -r] 
M*)= 2^i[3-a + e, -a, 2 - a, -r] 



where cq is determined by the condition V'( r )l r ^oc = ^- Note that r(— a) is singular when a = 1,2. 
We find that -0 cannot satisfy the condition V- , ( r )l r ^ 00 = when a — > 2, but there is no problem for 
a — ► 1. In this paper, we focus on the situation of a = 1 ~ 1.5 from the viewpoint of the N-body 
simulation. 




e=0.1 
a=1.05 
a=1.30 
ff=1.50 



0.0 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.S 1.0 

E/iAo E/(S 



Fig. 1. The energy part of the distribution function (left panel) and the differential energy distribution (right panel) 
of the New Generalized NFW model with the constant anisotropy parameter j3 = 1/2 and e = 0.1: dotted lines 
(o = 1.05, closes to the NFW model), dashed lines (a = 1.30), solid lines (a = 1.50, where the asymptotic behavior 
of density profile is similar to the Moore model). 



6 D. Ma and P. He 

4 p 1 1 1 1 1 1 




& 2 



EM, 



0.4 0.6 
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(1=1.5 
e=0.1 
e=0.6 

6=1.1 



Fig. 2. The energy part of the distribution function (left panel) and the differential energy distribution (right panel) 
of the New Generalized NFW model with the constant anisotropy parameter /3 = 1/2 and a = 1.5: dotted lines 
(e = 0.1, where the asymptotic behavior of density profile is similar to the Moore model), dashed lines (e = 0.6), solid 
lines (e = 1.1). 



Now we can get <?(?"#) from Eq. (|T2]) and Eq. (JTBJ): 



l{r E ) 



Bi(r) 



AitaT{~a)[B 2 {r)A A {r) + B 3 {r)A 5 (r)} 



Bx{r) = -(-2 + a - e) (-1 + a - e) [-1 + a + (2 + e)r] r (1 + r)~ 2+Q - £ 

B 2 (r) = (2 - a) + (4 + 2e - 3a)r + [(2 + e)(l + e) - a(3 + 2e)]r 2 

B 3 (r) = a - 2 + (2a - 3 - e)r 

^4(r) = 2 -Fi[l-a, 1-a + e, 2 - a, -r] 

A 5 (r)= 2-Fit-a, 1-a + e, 2 - a, -r] 



(17) 



Note that Bi(r) = if a — e = 1, 2. However, g{rE) doesn't approach zero when a — e — > 1,2. We 
also find that a should be larger than 1 as the distribution function should be positive. We couldn't 
get the Te from ip( r E) = E analytically but can get the numerical solutions. F(E,L) and dM/dE 
are expressed as: 



F(E,L) 



Bi(r) 



8tt 3 LaT(-a)[B 2 (r)A A {r) + B 3 (r)A 5 (r)} 



(18) 



dM _ r 2 Bi(r) 

~dE ~ 2aT(-a)[B 2 {r)A 4 {r) + B 3 {r)A 5 {r)} 



(19) 



Fig.[T]and Fig.[2]show g(E) — E/tJjq and dM/dE — E/ip for different values of parameters a and e. 
Here, i/)q — ip(r = 0). 
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3.2. Generalized Moore profiles 

The profile of the Moore modep] is shown as p oc l/((r/r ;s ) 1 - 5 (l + (r/rs) 15 )). We consider one 
generalized density profile of the Moore model: 

9 = ° (r/r s )«(l + (r/r s )^-<*) ' (20) 

As in section 3.1, we still set r s — 1, M tot — 1 and G = 1 here. We can not express C as 
a function of a and e analytically but can calculate C in a numerical way while considering the 
condition M tot = 4wr 2 p(r)dr = 1. 

The relative potential is 

</> (r) = C c p —7 -— 

(a — 2)(a — 3) 



3 — a + e 
d — a + e 

Ai(r) = 
- -3 + a 

3-^2 — ^— ) J-) 



(21) 



-2 + a 6-2a + e 5 - 2a + e 



„3-a+el 



-3 + a — e 3 — a 



where Co is determined by '0(r)| J ,_ >oo = numerically. tp{f) cannot satisfy the condition ■0( r )l r ^ O o = 
for some values of a > 2. Therefore, just as in section 3.1, we only focus on the situation of 
a = 1 ~ 1.5. Then g(r E ) can be calculated from Eq. (fT2|) and Eq. (|2Tj) : 



B 4 (r) 



Be(r) 



4 7 rS 5 (r)r 1 (-3 + a) A 7 {r) + T 2 A 6 {r) 



„3-a+ei 



(22) 



. , , ~ M -2 + a 5-2a + e 

A 7 {r) = 2 F 1 [1,—- ,- — , 

— 3 + a — e 3 — a + e 

S 4 (r) = (-l + a)r Q + (2 + e)r 3+£ 

B 8 (r)= (r Q +r 3+£ ) 2 

5 6 (r) = (6 - 5a + a 2 )r- 1+a 



In this kind of profile, we also calculate r E by numerical way. The F(E,L) and dM/dE are 
expressed as: 

B 4 (r) B 6 (r) 



F{E ' L) ^LB 5 {r)T 1 (-3 + a) A 7 (r)+T 2 A 6 (r) 



dM 
~dE 



r 2 B 4 (r) 



B e (r) 



(23) 



2B 5 (r)ri (-3 + a) A 7 {r) + T 2 A 6 (r) 
Fig. [3] and Fig. 0] show g(E) — E/ip and dM/dE — E/ip for different values of parameters a and e. 
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Fig. 3. The energy part of the distribution function (left panel) and the differential energy distribution (right panel) 
of the Generalized Moore model with the constant anisotropy parameter of /3 = 1/2 and e = 0.1: dotted lines (a = 1.05, 
where the asymptotic behavior of density profile is similar to the NF W profile's, but DF doesn't meet the condition 
dg(E)/dE > 0.), dashed lines (a = 1.30), solid lines (a = 1.50, closes to the Moore model). 
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Fig. 4. The energy part of the distribution function (left panel) and the differential energy distribution (right panel) 
of the Generalized Moore model with the constant anisotropy parameter of [3 = 1/2 and a = 1.3: dotted lines (e = 0.1), 
dashed lines (e = 0.6), solid lines (e = 1.1, where DF doesn't meet the condition dg(E)/dE > 0). 



Table 1. Parameters e and a c 



c = 


0.01 


0.1 


0.2 


0.3 


0.4 0.5 


Ctc ~ 


1.110 


1.127 


1.146 


1.166 


1.186 1.209 


c = 


0.6 


0.7 


0.8 


0.9 


1.0 


Ctc ^ 


1.230 


1.253 


1.276 


1.299 


1.323 



In this model, a basic stability conditiorPSESETl f or spherical stellar system dF/ dE > is 
equivalent to the condition dg(E)/dE > (note: the " E" here is not the energy of a particle but the 
binding energy). We note that spherical systems with certain a and e may not meet this condition. 
This condition requires that the parameter should be larger than a critical value a c for the given 
parameter ej. We calculate a c numerically and show them in Table [T] The difference of the two 
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generalized profiles is shown in Fig. [5l One can see that the parameter a controls the asymptotic 
density profile in the inner parts, while e is responsible for the asymptotic behavior of the outer 
profile. 




Fig. 5. We set G = 1, Mtot = 1 and /3 = 1/2. Two solid lines: New Generalized NFW profiles with different 
parameters. Two dashed lines: Generalized Moore profiles with different parameters. 



4. Discussion and Conclusions 

N-body simulations show that the anisotropy parameter (3 « 0.5 in the outer part of dark matter 
halos. The distribution functions with constant anisotropy parameter /3 = 1/2 and the corresponding 
differential energy distributions can be calculated from the formulae (fT3i|) . N-body simulations also 
show that the logarithmic slope of density profile j = 1 <~ 1.5 in the inner part of the dark matter 
halo (p ~ r~ 7 ). 

Many authors work at construction of the dark matter halo models which are better to be 
analytical, simple, realistic and having flexible anisotropy profile. Unfortunately, it is difficult to 
meet these requirements at the same time. 

In this paper, we only consider the models with constant anisotropy parameter (3 — 1/2 and 
then calculate the DFs and DEDs of the two models, the new generalized NFW model and the 
generalized Moore model. Both have two parameters e and a for determination of the asymptotic 
behavior of density profile in the inner and outer parts of the halo. With these two parameters, our 
models can cover many relatively realistic density profiles which come from N-body simulations. A 
physical postulate requires that DFs should be positive in the phase space, and our DFs satisfy this 
basic condition. A basic stability condition for a spherical stellar system dF/dE > is the sufficient 
condition for the isotropic case, but is no more a sufficient condition for the anisotropic onePSEHEZl 
However, we still can use this condition to do some rough discussion. In the new generalized NFW 
model, all DFs with e > and a — 1 ~ 1.5 meet the condition dF/dE > 0. In the generalized Moore 
model, we find that not all DFs with e > and a = 1 ~ 1.5 satisfy this condition. In order to satisfy 
the condition dF/dE > 0, on should be larger than a critical value a c for a given parameter e;. 
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We have not dealt with the stability of our models in details in our paper. Although it is difficult 
to determine the stability domain in the parameter space of the halo model, the study of stability 
of th e stellar sys tem and dark matter halo is very important and some works have already been 
done I 38 | 39 | 4(l | 41 | To be 

more realistic, a more complicated case, the axisymmetric model should 
be considered^ 43 ' 44 ' 45 '. Construction of the polycomponent models is also necessary as the stellar 
component or the central black hol e ' ' ' usually combines with the dark matter component 
and t he two-component models with the stellar and dark matter components were constructed by 
Ciottffl 

Although the density profiles of our models are realistic enough and we assume that the anisotropy 
parameter /? = 1/2 from the vie wpoint of the simulation, it is still not realistic enough for the 
anisotropy profile. An & Evans^ explored the model whose DF has the form as below, 

F(E,L)=Y,L- 2 ^ME). (24) 

i 

DF is superposition of two or more terms here and this case is more complicated, nevertheless, this 
kind of DF has more flexible anisotropy profile and it is revealing for our further work to pursue a 
greater variety of anisotropic behavior of the dark matter halo models. 
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Appendix A. Hypergeometric Function^ 

pF q is the regularized hypergeometric function which is defined as 

TTi/ l t \ pFq( a i> ■■■! a p] b\, b q ; z) , . 

p F q (ai, ...,a p ;6i, ...,b q ;z) = r ^ r ^ ^ , (A.l) 

where T(z) is a gamma function. And p F q (ai, a p ; &i, b q ] z) is the generalized hypergeometric 
function: 

p F q { ai , a p ; b u b q - z) = ^ 77 s 77^ ttt- tt , (A.2) 

^ {0\)k (02)k •■■ \0 q )k «! 

where (a)k is the Pochhammer symbol, 

, . r(a + k) , , . s 

a ) fc = v , ' =aa+l ... a + fc-1 . A.3 
T(a) 

The specific hypergeometric functions can be calculated by mathematical software packages such as 
Mathematica. 



a See http: //mathworld. wolfram. com/ 
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